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Abstract 

The explicit criterion is derived in detail for the convergence of the Magnus 
expansion and the existence of the Magnus solution in the interaction picture, 
i.e., the exponential propagator in the weakly coupled spin (1=1/2) system 
S„AMX... (n=l,2,3,...) in which only spin group S is subjected to a time- 
dependent shaped selective radiofrequency pulse. The derivation is built on the 
same scheme as the Maricq's (J. Chem. Phys. 86 (1987) 5647). It is shown that 
the criterion depends only upon amplitude of the time-dependent field applied 
to the system, and the Magnus expansion converges and the Magnus solution 
exists when the flip angle of a non-negative-amplitude shaped RF pulse or a 
weak-amplitude shaped pulse is smaller than 2tt. The exponential propagator 
then can be decomposed into a product of a series of elementary propagators and 
can be used to determine time evolution of the spin system under the shaped 
pulse. The linear differential equations to determine the unknown parameters 
in the Magnus solution are obtained explicitly. An alternative propagator in an 
expansion form for the Magnus solution is also proposed to describe the time 
evolution when the criterion is not met. 
PACS: 03.65.Ge; 33.25. +10 
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The Magnus expansion [1] plays an important role in determination for time 
evolution of quantum systems under time-dependent Hamiltonian. The advan- 
tage of the Magnus expansion over the standard quantum mechanical perturba- 
tion theory is that Hermiticity is retained to each order in the expansion. The 
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Magnus expansion has an extensive application in nuclear magnetic resonance 
(NMR) spectroscopy: the average Hamiltonian theory [2, 3] which is based on 
the Magnus expansion has been used extensively in the design and analysis of a 
variety of multipulse NMR experiments both in liquid and in solid states. The 
Magnus expansion has also been used successfully to investigate the response 
of coupled spin systems to a time-dependent radiofrequency (RF) field, e.g., 
a shaped selective pulse [4, 5]. The Magnus solution, namely, the continuous 
exponential solution to the time-dependent Schrodinger equation, has recently 
been used to characterize the behavior of a coupled spin system under the time- 
dependent radiation damping field [6]. Although the Magnus expansion and 
the Magnus solution are applied conveniently to a variety of time-dependent 
systems, there are strict limits on the existence of the solution and there is also 
a possible serious problem for the expansion concerning its convergence [1, 7 - 
9]. As pointed out by Magnus [1] as long ago as in 1954, the Magnus solution 
U(t) = exp[— iil(t)} exists only when the so-called Magnus criterion is satisfied: 
|A,(t) — Xj(t)\ 7^ n2n, where Aj(i) are the eigenvalues of the operator fi(f) and n 
is a nonzero integer. This criterion is an implicit one and has little value since 
Aj(t) are the unknown quantity under investigation. An explicit criterioi^ is 
therefore needed for the convenient application of the Magnus solution. Maricq 
[8] showed that the operator fi(t) exists for short evolutional times and the Mag- 
nus expansion always converges in the neighborhood of the time t = 0. Maricq 
[9] further found that the convergence for an isolated two-level system under a 
linear polarized field can be improved: the Magnus expansion is convergent and 
hence the Magnus solution exists in time interval [0, r] if the following inequality 
is satisfied 



where ±w(t)/2 are the two eigenvalues of the time-dependent Hamiltonian H(t) 
of the system. The advantage for this criterion is that it is dependent only on 
the eigenvalues of the Hamiltonian H(t). However, an explicit criterion may 
be far from the simple form of inequality (1) and has not yet been found for a 
general time-dependent multilevel system [7 - 9]. 

A shaped selective pulse is an amplitude- modulation and/or phase- modulat- 
ion RF field [4, 10 - 13]. The Magnus solution may provide a convenient starting 
point for an exact description of time evolution of coupled spin systems under a 
time-dependent shaped selective RF pulse. However, so far an explicit criterioi^ 
for existence of the Magnus solution and convergence of the Magnus expansion 
has not 6eer0 established for the general spin systems. This general problem 
could not be settled simply, but in some special cases it may be able to be solved 
without much difficulty. In this communication the Maricq's criterion (1) for a 
single two-level system is extended to the weakly coupled spin (1=1/2) systems 

3 The original words are modified here. The original word(s) appears in the original 
manuscript and its published paper (Xijia Miao, Phys. Lett. A 271 (2000) 296-302). Hereafter 
use italics to mark the modified word(s). 
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S„AMX... (n=l,2,...), where only spin S is undergone a time-dependent RF 
field. It is known that the properties of the Magnus expansion depend strongly 
on the representation [7 - 9]. The present criterion is derived explicitly in the 
interaction picture defined by spin Hamiltonian of the spin system in the absence 
of the time-dependent RF held. In this representation the criterion is dependent 
only upon the amplitude of the time-dependent RF held applied to the system. 
Therefore, in practice this is quite convenient to use the criterion. 

For simplification, consider the SAMX... weakly coupled spin (1=1/2) sys- 
tem in which only spin S is subjected to a time-dependent RF field, e.g., a 
shaped selective pulse. The contribution of the time-dependent RF field to spin 
Hamiltonian of the system is given by 

Hi(t) = uJi(t){S x cos <j>{t) + S y sin <j>{t)} (2) 

where 0J\(t) and <p(t) are known as the amplitude and phase of the time- 
dependent RF field, respectively. Then the total spin Hamiltonian of the spin 
system in a strong static magnetic field can be written in the rotating frame 

H{t) = H a + H 1 {t) (3) 

where H = Hj + CliS z \ Hi = J^k &khz + J2k<i ^Jkl^hzhz and tti = Q s + 
2 J2k ^Jkshz', the symbol I denotes the spin A, M, X, .... Obviously, both the|j 
operators Hi and fli are members of the longitudinal magnetization and spin 
order (LOMSO) operator subspace of the Liouville operator space of the spin 
subsystem AMX... of the system SAMX... [14]. They Q commute with each 
other because any two base operators of the LOMSO operator subspace com- 
mute with each other [14] . They also commute with the magnetization operators 
{S Xl S v , S z } of spin S. The time-evolutional propagator corresponding to the 
time-dependent Hamiltonian (3) can be written as 

U(t) = Texp(-i / H{t')dt') = exp(-iH t)Ui(t) (4) 
Jo 

where T is Dyson time-ordering operator and the propagator Ui(t) in the inter- 
action representation is given by 

Ui{t) = TcxpH / Hi{t')d£) (5) 
Jo 

with the time-dependent Hamiltonian in the interaction representation defined 
by the unperturbed Hamiltonian Hq of Eq. (3): 

Hi(t) = expiiHo^H^t) exp(-iH t) (6) 

By expanding the right-hand side of equation (6) according to the rotation 
transformations [2] between any two product operators and the fact that [Hi, 
Qi] = 0, the Hamiltonian Hi(t) is reduced to an explicit expression 

H^t) = un(t)h(t).S = wi(t){ftx(t)Sx + h y (t)S y + h z (t)S z } (7) 

6 Here drops a redundant original word. 
7 Here drops a redundant original word. 
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where the operator h x (t) = cos(— fljt + 4>{t)), h y (t) = sin(— Ujt + <f>(t)), and 
h z (t) = 0. Because h x (t) 2 + h y (t) 2 + h z (t) 2 = E (unit operator), the opera- 
tor h(t) = (h x (t),h y (t),h z (t)) is an operator vector with the unit magnitude. 
Obviously, the operator h(t) belongs the LOMSO operator subspace of the 
subsystem AMX.... It follows from Eq. (7) that 

\mt)\\<\^(t)\-\\Ht)\\.\\s\\ (8) 

Therefore, all eigenvalues of the Hamiltonian -ff»(i) are bounded by the am- 
plitude in magnitude of the time-dependent RF field. In a single-spin 
system h(t) is a vector of simple scalars and in this case the criterion, that is, in- 
equality (1), was obtained explicitly by Maricq [9], but in the SAMX... coupled 
spin system h(t) can be taken as a vector of the operators h p {t) (p = x,y,z) 
that all commutes one with each other because the operators h p (t) are members 
of the LOMSO operator subspace. In the present more general case the crite- 
rion of the Magnus expansion in the interaction picture has not been obtained 
explicitly and it is just the purpose of the present work to extend the criterion 
(1) to the S„AMX... coupled spin systems. 

In order to determine exactly and analytically time evolution of the cou- 
pled spin system under a time-dependent RF field one needs to solve the time- 
dependent Schrodinger equation or equivalent Liouville-von Neumann equation. 
However, it is usually difficult to solve exactly the time-dependent dynamical 
equations and therefore to determine exactly the density operator and the re- 
sponse of the system to the time-dependent RF field. A numerical calculation is 
a general routine to solve time-dependent problems, but it is usually not much 
helpful for one to gain a deeper insight into the physical essence for the effect 
of a time-dependent RF field on the time evolution and for the effect of scalar 
coupling on the response. The approximated treatments that usually employ 
the perturbation theory or the average Hamiltonian theory are valuable, but 
there are some limitations and disadvantages of these approximated methods 
and therefore, an exact solution to the time-dependent problems is still attrac- 
tive. One strategy to determine exactly and analytically time evolution of the 
coupled spin system under a time-dependent RF field is that the propagator 
Ui(t) of Eq. (5) is written as an exponential operator without the Dyson time- 
ordering operator T, as suggested by Magnus in early day [I]. Note that the 
magnetization operators {E, S x , S y , S z } of the spin S form the u(2) Lie algebra 
subspace and the longitudinal magnetization and spin order operators {E, Ikz, 
2Ikzhz, ^Ikzhzlmz, ■■■} of the subsystem AMX... form the LOMSO operator 
subspace [14] by including the unit operator E. Then the complete operator set 
{E, S p , Ikz, 2S p Ikz, 21kzllz, ^Splkz Iiz, p = x,y, z} is a closed operator al- 
gebra subspace of the Liouville operator space of the spin system SAMX... and 
actually, this subspace is the direct product space of the subspaces u(2) and 
LOMSO, that is, G = u{2) LOMSO [14]. Obviously, the Hamiltonian H^t) 
of Eq. (7) is a member of the operator subspace G because the operators h p (t) 
(p = x,y,z) in Eq. (7) are members of the LOMSO operator subspace. Then 
the propagator Ui(t) of Eq. (5) corresponding to the time-dependent Hamilto- 
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nian Hi{t) is also a member of the subspace G [14] and may be conveniently 
written as an exponential unitary operator: 

Ui{t) = exp[-ifi(i)] (9) 

with the Hermitian operator 

n(f) - n x (t)s x + Q y (t)s y + n z (t)s z (10) 

where the Hermitian operator tt(t) is an element of the operator subspace G 
and Q p (t) (p — x, y, z) belong the LOMSO operator subspace. The exponential 
propagator of Eq. (9) is quite simple and can be further decomposed into a 
sequence of elementary propagators [15, 16] 

Ui(t) = exp(-ia(t)S z ) exp(-0(t)S y ) exp(-iCl(t)S z ) 

x exp{i/3(t)S y ) exp(ia(t)S z ) (11) 

where the operators a(t), f3(t), and Q(t) that are also members of the LOMSO 
operator subspace are related to the operators Cl x (t), Q y {t) 7 and Q z (t) in Eq. 
(10) by the following formulae: 

tl x (t) sin a(t) = Q y (t) cos a(t) (12a) 

n z (t)sm/3(t) = [tt x (t) cos a(t) +Q y (t) sin a(t)] cos 0(t) (12b) 

Cl(t) = tt z (t) cos f}{t) + [Cl x (t) cos a(t) +D, y (t) sin a(t)} sin (3 (t) (12c) 

The exponential propagator of Eq. (9) of the SAMX... spin system can be 
extended to the coupled spin systems S n AMX... (n=l,2,3,...) and in this more 
general case the propagator can be written at@ 

n 

Ui(t) - I] + toy(t)Sk» + ttz(t)S kz )} (13) 

fc=i 

With the help of the decomposed propagator of Eq. (11) one can determine 
conveniently time evolution of the coupled spin system SAMX... under a shaped 
pulse. However, there are two problems that need to be solved in advance. One 
of which is that the existence of the exponential propagator Uiif) of Eq. (9) 
needs to be verified. This would be achieved by setting up the convergence 
condition of the Magnus expansion of the operator f2(t) in Eq. (9) and hence 

8 The incorrect expression (13) appearing in the original manuscript and its published paper 

is 

n 

Ui(t) = n exp[-in x (t)s kx + n y (t)s ky + n z (t)s kz ] (13) 
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the existence condition of the Magnus solution of Eq. (9) in the interaction 
picture. Another is that the unknown parameters in the decomposed propagator 
of Eq. (11) should be determined explicitly. It was shown by Maricq [9] that the 
explicit condition (1) for the convergence of the Magnus expansion in a single- 
spin system under a linear polarized field is related only to the eigenvalues of 
the time-dependent Hamiltonian of the system. However, it has not been clear 
whether the convergence of the Magnus expansion in the coupled spin s^errf! 
SAMX... is still dependent only upon the eigenvalues of the time-dependent 
Hamiltonian of the system. From the point of view of convenient application, 
the convergence had better depend only on the eigenvalues of the Hamiltonian 
or even only upon the amplitude and/or phase of the time-dependent RF field 
applied to the system in the interaction frame. On the other hand, it is expected 
that the convergence of the Magnus expansion in the interaction frame may 
depend only upon the amplitude of the time-dependent RF field as all the 
eigenvalues of the time-dependent Hamiltonian (6) in the interaction frame are 
bounded on by the amplitude of the RF field in magnitude (see inequality (8)). 
Therefore, in the following the convergence criterion is derived in the interaction 
representation enabling it to depend only upon the amplitude of the RF field. 
The derivation is based on the decomposed propagator of Eq. (11) and the 
time-dependent Schrodinger equation in the interaction frame: 



The decomposed propagator of Eq. (11) is first expanded as a linear combination 
of the magnetization operators S p (p — x, y, z) 



Equation (15) is the expansion expression of the propagator Uj(t) in the operator 
subspace G. Then inserting both the Hamiltonian Hi(t) of Eq. (7) and the 
expansion-form propagator Ui(t) of Eq. (15) into equation (14) one obtains a 
set of equations entirely equivalent to the time-dependent Schrodinger equation 




iHi(t)Uj(t) (h=l) 



(14) 




(15) 




= -wi(t)cos(-ft(i))[fc x (i), h y {t), h z (t)} 
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+-Wx{t) sm(-n(t))[-h z (t) sin a(t) sin (3(t) + h y (t) cos (3(t), 
h z (t) cos a(t) sin/3(t) — h x (t) cos j3(t), 

h x (t) sina(i) sin f3(t) — h y (t) cos a(t) sin (3 (t)] (16b) 

Equations (160 are a set of linear coupled differential equations with the opera- 
tor variables: cos(|f2(i)), cosa(i) sin/3(i) sin(|fi(t)), sin a(t) sin j3(t) sm(hCl(t))X^\ 
and cos/3(t) sin(|0(t)) of the LOMSO operator subspace. They can be used to 
determine generally the unknown parameters in the decomposed propagator of 
Eq. (11). Also, from them one may extract the criterion of convergence of the 
Magnus expansion in the interaction frame. Because any two base operators 
of the LOMSO subspace commute and the operator Cl(t) belongs the LOMSO 
subspace one has 

| cos{\m) = -i sm(ifi(t))|fi(t) (17) 
Noting equation (17), equation (16a) leaves one a simple form 

^-Cl(t) =wt(t){h x (t)cosa(t) sxn0(t) 



+h y (t)sma(t)smP(t) + h z (t) cos/?(t)} (18) 

The equation (18) can be integrated conveniently because all the operators in 
the equation belong the LOMSO subspace. By employing the initial condition 



10 The incorrect equations (16) appearing in the original manuscript and its published paper 
are 

— cos(-n(t)) = -oii(t) sin(-£l(t)){h x (t) cos a(t) Bin 0(t) 

+h y (t) sin a(t) sin 0{t) + h z (t) cos 0(t)} (16a) 

— [cos a(t) sin /?(t) sia(if2(t)), sina(t) sin^(t) sin(if2(t)), cos^(i) sin(ifi(t))] 
dt 2 2 2 

= -u>i(t)cos(if2(t))[M*), h y (t), h z (t)} 

+-u)l (t) sin(-f2(t)) [h z (t) sin a(t) sin 0{t) - h y (t) cos 0(t), 
h z (t)cosa(t) sin/3(t) ~ h x (t) cos/3(t), 

hx (t) sina(t) sin 0(t) - h y (t) cos a(t) sin 0(t)] (16b) 

llr The incorrect operator variable appearing only in the published paper is sino(t) sin/3(i) : 
sin(if2(i)) 
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Cl(0) = the integration is written at0 

Cl(t) = [ di>i(t')M*')-9(*')} (19) 
Jo 

where the operator vector with unit magnitude g(t) = (cosa(i) sin f3(t), sina(i) 
x sin/3(i), cos/3(i)). 

Now the eigenvalues of the operator fi(i) in the Magnus solution (9) are 
readily calculated in the direct product basis of eigenbase \i) and s) of single- 
spin angular momentum operators Ikz and S z by diagonalizing the operator 
Q(t) (see Eq. (11)) 

Xt(t) = (i\(* a \Cl(t)\* a )\i) = (i\(a\fl(t)S,\a)\i) = m s (i\fl(t)\i) (20) 

where the magnetic quantum number m s = +1/2 and —1/2 for s = and 
1, respectively. Because of the unit magnitude of the operator vectors h(t) 
and g(t) it follows from Eq. (19) that each eigenvalue of the operator Cl(t) is 
bounded by 

|(i|fi(t)|i)| < / dt'\wi(t')\ (21) 

This inequality (21) together with equation (20) can lead to directly the desired 
result by the triangle inequality: 

\K{t)-^{t)\ < \\m\ + \^(t)\ < f dt'\u x {t')\ (22) 

Jo 

By comparing the inequality (22) with the Magnus criterion one obtains the 
explicit criterion of convergence of the Magnus expansion in the interaction 
frame: 

I(t) = / diViCOl < 27T ( 23 ) 
Jo 

The criterion (23) is the same as the Maricq's criterion (1) and its derivation is 
also built on the same scheme as the Maricq's in the case of single-spin systems. 
However, it should be pointed out that in the Maricq's derivation both h(t) and 
g(t) in Eq. (19) are the vectors of simple scalars, while in the present situation 
both the two vectors h(t) and g(t) are the operator vectors of the LOMSO 
operator subspace. 

When a time-dependent shaped RE pulse such as Gaussian pulse [10] and 
sech pulse [6] applied to the coupled spin system has a positive amplitude 
wi(t) > over the whole pulse duration [0, t] the integration I(t) equals ac- 
tually flip angle of the shaped RF pulse. Then the Magnus expansion in the 
interaction frame converges and the Magnus solution of Eq. (9) exists when 

12 The incorrect equation (19) appearing only in the published paper is 

Cl(t)= P dt'{uj 1 }(t')h(t').g(t') (19) 
Jo 
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flip angle of such shaped pulse is smaller than 27r, that is, the criterion (23) 
is met. Evidently, the explicit criterion (23) serves as a general lower limit for 
the convergence radius of the Magnus expansion in the interaction picture [9] . 
Thus, that the condition (23) is not fulfilled does not necessarily imply that the 
Magnus expansion diverges and the Magnus solution of Eq. (9) does not exist. 
The Magnus solution exists really as long as the Magnus criterion is satisfied 
[1]. However, the Magnus expansion always converges and hence the Magnus 
solution, i.e., the exponential propagator Ui(t) of Eq. (9) in the interaction pic- 
ture exists certainly when the criterion (23) is met. Therefore, the exponential 
propagator of Eq. (9) and the decomposed propagator of Eq. (11) can be used 
to determine exactly time evolution of the coupled spin system S„AMX... under 
a non-negative-amplitude shaped RF pulse such as Gaussian [10] and sech [6] 
pulses applied to the spin S when flip angle of the shaped RF pulse is smaller 
than 2n, a very popular situation in high-resolution nuclear magnetic resonance 
spectroscopy. 

For a weak-amplitude time-dependent RF field the dot product of the two 
operator vectors h(t) and g{t) in equation (19) is usually nearly equal to unit 
operator. As a consequence, each eigenvalue of the operator £l(t) is evaluated 
approximatively b}Q 

(i\Cl(t)\i) » / dt'wi(i') (24) 
Jo 

Then the explicit criterion (23) is reduced to the following inequality 

9(t) = [ dt'uj^t') < 2tt (25) 
Jo 

13 Here drops a redundant original word. 

14 The approximate formula (24) is not easy to understand. When the amplitude wi(t) of 
an RF field is small, the norm that is proportional to |u>x(i)| is also small, as can be 

seen in (8). Then in that case the Magnus expansion is quickly convergent [8] and thus one 
may use reasonably the lowest-order approximation formula to calculate the operator Q(t), 

Q(t) ~ f dt'Hi(t'). 
Jo 

Now the operator tt(t) may be diagonalized by using the unitary transformation obtained 
from Eqs. (9) and (11). On the other hand, the interaction Hamiltonian Hi(t) also can 
be diagonalized by using the unitary transformation obtained from Eqs. (2) and (6). One 
therefore obtains from the lowest-order approximation formula that 

Cl{t)S z ra f dt'{ui(t')V(t,t') + S z V(t,t')}. 
Jo 

Here V(t, t') is a unitary operator. In the trivial case that the operator V(t, t') is the unit 
operator the formula (24) holds naturally. In a general case the unitary operator V(t, t') may 
be written as 

V(t,f) = exp{-i(V^(t,t')S x + V*(t,t')Sy + V/(t,i')S z )}. 

Here Vj(t, t') (p = x, y, z) are the LOMSO operators of the subsystem AMX.... Note that the 
amplitude function D\(t) and phase function <f>(t) in (2) can be taken as arbitrary functions, 
respectively. Then the lowest-order approximation formula shows that one must set V/(t, t') « 
and Vy (t, t') ss 0. This directly leads the lowest-order approximation formula to the formula 
(24). 



9 



This criterion shows that the Magnus expansion converges and the Magnus 
solution (9) exists when flip angle 9(t) of a weak-amplitude time-dependent RF 
pulse is smaller than 2ir, that is, the inequality (25) is met. The criterion (25) 
may be useful for the case of those shaped RF pulses such as sine and hermite 
pulses [4] whose amplitude u>i(i) may not be always positive over the whole 
pulse duration. 

The real unknown parameters in the propagators of Eqs. (9), (11), and (15) 
can be generally determined by solving the linear differential equations (16). It 
should be noted that the linear differential equations (16) are not independent 
but in which the four operator variables are constrained 15 ] by 

[cos(ifi(t))] 2 + [cos a(t) sin /3(f) sin(ift(t))] 2 

+[sina(fjsm/3(t) sin(-fi(f))] 2 + [cos/?(t) sin(-n(t))] 2 = E 

where E is unit operator. 

Actually, the propagator Ui(t) of Eq. (5) can be generally expressed as a 
sum of the base operators of the operator subspace G 

U I (t) = f(t)-2iJ29 P (t)S P (P = x,y,z) (26) 
p 

where the operators f(t) and g p {t) are the LOMSO operators. Because of the 
unitarity of the propagator Ui(t) the operators are written respectively in form 
as fit) = cos(|fi(t)); g x (t) = cosa(i) sin/3(i) sin(|f2(t)), g v (t) = sina(t)sin/3(t) 
x sin(|0(i)), and g z {t) = cosj9(i)siri(if2(t)), where a{t), f3(t), and tl(t) are 
still the LOMSO operators. Obviously, a(t) = a(t), /3(t) = /3(t), and U(t) = 
Ct(t) when the Magnus solution (9) exists. By inserting the expansion-form 
propagator of Eq. (26) into the Schrodinger equation (14) one can still obtain the 
same linear differential equations as equations (16) but with the new operator 
variables f(t) and g p (t) (p — x,y,z). The differential equations are entirely 
equivalent to the Schrodinger equation (14) and can determine an alternative 
solution for the Magnus solution (9) to the time-dependent Schrodinger equation 
[9]. Because the solution to the differential equations always exists one can 
obtain the unknown parameters in the propagator of Eq. (26) by solving the 
differential equations. Thus, the expansion-form propagator of Eq. (26) always 
exists and can be generally used to determine time evolution of the coupled 
spin system under an arbitrary time-dependent RF pulse no matter whether 
the criterion (23) is met or not. This indicates that the time evolution can be 
described in the small operator algebra subspace G = u{2) (^) LOMSO instead 

of the whole Liouville operator space of the system. Thus, this simplifies 
greatly determination of the time evolution. 

15 The original manuscript's word is corrected here. 
16 Here drops a redundant original word. 
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In general, the integration I(t) is always more than flip angle 9(t) for a 
shaped pulse. Then before employing the criterion (23) one usually needs to 
first calculate explicitly the integration I(t) for complex shaped selective pulses 
such as the BURP pulses [11] and the Gaussian pulse cascades [12, 13] whose 
amplitude functions uj\(t) are not always non-negative over the whole pulse 
duration even if their flip angles given in advance are smaller than 2tt. The 
numerical calculation^! shows that the criterion (23) is met for the excitation 
pulses G 4 and Q 5 of the Gaussian pulse cascades [12, 13] but is violated for the 
BURP pulses [11] including E-BURP, U-BURP, I-BURP, RE-BURP pulses and 
the Gaussian pulse cascades G 3 and Q 3 [12, 13]. It should be pointed out that 
the violation of the criterion (23) for these pure-phase shaped selective pulses 
does not imply that the exponential propagator of Eq. (9) can not be used to 
describe time evolution of the coupled spin system under these shaped pulses. 
However, as an alternative one still can use the expansion-form propagator of 
Eq. (26) of the smaller operator subspace G = u(2) LOMSO instead of that 
one of the whole Liouville operator space of the coupled spin system S n AMX... 
to determine generally time evolution of the spin system under those shaped 
pulses and shaped decoupling sequences [17] that do not satisfy the criterion 
(23). 
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